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ABSTRACT 

Grover’s  oracle  based  unstructured  search  algorithm  is  often  stated  as  “given  a  phone  number  in  a  directory,  find  the 
associated  name.”  More  formally,  the  problem  can  be  stated  as  “given  as  input  a  unitary  black  box  Uf  for  computing  an 
unknown  function /{ 0,1}"  — >{0,1 } find  x=x0  an  element  of  {0,1}''  such  that  f[x0)  =1,  (and  zero  otherwise).  The  crucial 
role  of  the  externally  supplied  oracle  Uf  (whose  inner  workings  are  unknown  to  the  user)  is  to  change  the  sign  of  the 
solution  |x0)  ,  while  leaving  all  other  states  unaltered.  Thus,  [/depends  on  the  desired  solution  x0.  This  paper  examines 

an  amplitude  amplification  algorithm  in  which  the  user  encodes  the  directory  (e.g.  names  and  telephone  numbers)  into  an 
entangled  database  state,  which  at  a  later  time  can  be  queried  on  one  supplied  component  entry  (e.g.  a  given  phone 
number  t0 )  to  find  the  other  associated  unknown  component  (e.g.  name  Xo).  For  N=2"  names  |  x)  with  N  associated  phone 

numbers  1 1) ,  performing  amplitude  amplification  on  a  subspace  of  size  N  of  the  total  space  of  size  N2  produces  the 
desired  state  |x0)|t0)  in  VlV  steps.  We  discuss  how  and  why  sequential  (though  not  concurrent  parallel)  searches  can  be 
performed  on  multiple  database  states.  Finally,  we  show  how  this  procedure  can  be  generalized  to  databases  with  more 
than  two  correlated  lists  (e.g.  |x)|t)|,s)|r) ...). 

Keywords:  quantum  computing,  Grover’s  search  algorithm,  quantum  algorithms 


1,  INTRODUCTION 

In  addition  to  Shor’s  factorization  algorithm1,  Grover’s  search  algorithm2  is  another  highly  recognized  quantum 
algorithm,  being  widely  taught  in  many  texts  on  quantum  computation,3’4  and  serves  as  a  benchmark  for  nascent  physical 
implementations  of  quantum  computers5.  Grover’s  search  algorithm  (GSA)  considers  the  following  scenario6,  suppose 
you  have  a  large  table  T[0. .  ./V- 1 1  of  N  entries  for  which  you  would  like  to  find  some  element  z0.  More  precisely,  you 
wish  to  find  an  integer  x0  such  that  0  <  x0  <  N  and  T[x0]=  z0,  provided  that  such  an  x0  exists.  If  the  table  is  sorted  the 
problem  can  be  solved  in  a  time  Oflog  N).  Flowever,  in  many  interesting  problems,  ordering  or  structuring  the  data  may 
not  be  possible  or  practical,  and  one  must  resort  to  the  brute  force  method  of  exhaustively  searching  through  all  the  data 
until  the  result  is  found  (or  to  determine  if  it  even  exists).  Classically,  there  is  no  algorithm  that  succeeds  with 
probability  greater  than  !4  without  searching  through  more  than  half  the  entries  of  T.  Grover1  described  his  algorithm  as 
finding  a  needle  in  the  haystack,  and  equivalently  as  finding  the  associated  name  in  a  telephone  book  when  one  is 
supplied  with  a  given  telephone  number  (in  which  the  telephone  book  is  sorted  on  the  names,  but  random  on  the 
telephone  numbers).  Grover’s  quantum  unstructured  search  algorithm  can  solve  this  problem  on  a  quantum  computer  in 
expected  time  in  0(^N).  The  GSA  has  also  been  shown  to  be  optimal6,  implying  that  a  quantum  algorithm  cannot 
achieve  faster  than  a  quadratic  speedup  over  its  classical  counterpart. 

The  GSA  utilizes  an  oracle,  which  computes  a  function  /(x)  of  the  input  x,  but  whose  inner  workings  is  unknown  and 
unavailable  to  the  user.  The  Grover  search  problem  can  be  stated  formally3  as 

The  Grover  Search  Problem 

Input:  A  black  box  (oracle)  Uf  for  computing  an  unknown  function/:  {0,1}''  —>{0,1}. 

Problem:  Find  an  input  x0  {0,1}''  such  that/(x0)  =1  and  /(x/x0)  =0. 

In  the  above, /is  the  classical  function  which  evaluates  to  “yes”  on  the  needle  and  “no”  on  the  more  abundant  pieces  of 
hay  in  the  haystack.  Uf  is  the  unitary  representation  of  f  which  acting  on  x  encoded  into  a  quantum  multi-qubit  state  |x) , 
performs  the  reversible  operation 
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Uf\x)®\y)  =  \x)®\y®f(x)),  (1) 

where  |  y)  is  a  single  auxiliary  qubit  and  ®  denotes  binary  (mod  2)  addition.  (Note:  from  now  on  we  will  often  write  the 
tensor  products  of  state  |  x)  ®  |  y'j  simply  as  |  *)  \  y)  =  \  x,  y)  )• 

As  is  well  known,  and  as  will  be  explicitly  illustrated  below,  Uf  requires  knowledge  of  the  solution  x0  in  order  to  be 
explicitly  constructed4.  This  is  why  the  oracle  Uf  is  part  of  the  input  to  the  GSA,  and  it  has  to  be  supplied  externally  to 
the  user  performing  the  search.  Recently,  there  has  been  interest  in  developing  algorithms  that  would  dispense  with  the 
Grover  oracle  [7^  and  encode  the  search  list  directly  into  a  quantum  database  state  which  can  be  initially  constructed  (e.g. 
an  encoding  of  a  telephone  book),  and  subsequently  searched  at  a  later  time  (e.g.  given  a  telephone  number,  find  the 
associated  name).  Xu  et  al.1  have  designed  such  an  0(^lN)  algorithm  based  on  adiabatic  quantum  computing  (AdQC)  and 
experimentally  demonstrated  its  operation  on  a  two  qubit  ‘‘telephone  book”  in  an  NMR  quantum  computer.  In  their 
work,  only  the  names  were  encoded  into  the  quantum  database  state,  while  the  telephone  numbers  were  encoded  as 
classical  integers.  The  goal  of  this  paper  is  enunciate  a  quantum  search  algorithm  (QSA),  analogous  in  spirit  to  Xu  et 
al.1 ,  but  in  the  usual  quantum  circuit  model  paradigm  (i.e.  an  explicit  unitary  operator  approach  vs  the  Hamiltonian 
approach  of  AdQC). 

In  the  following,  we  briefly  recall  the  GSA  in  a  now  standard  form  that  elucidates  its  functionality  and  its  generalization 
to  quantum  amplitude  amplification  (QAA)  algorithms.  In  addition,  we  illustrate  the  GSA  on  a  simple  two  qubit  search, 
explicitly  constructing  the  Grover  oracle  Uf  revealing  its  dependence  on  the  solution  xq.  We  then  briefly  describe  the 
non-oracle  AdQC  database  search  of  Xu  et  al.1  and  describe  why  it  succeeds.  We  illustrate  their  algorithm  (as  their  paper 
does)  on  a  simple  telephone  book  search  example.  Next,  we  turn  to  the  main  portion  of  this  paper,  in  which  we  detail  a 
QSA  where  both  the  names  and  telephone  numbers  are  encoded  into  a  quantum  database  state.  For  N=2"  names  with 

N  associated  phone  numbers  |  /)  we  show  how  our  QSA  performs  amplitude  amplification  on  a  subspace  of  size  N  of  the 
total  Hilbert  space  of  size  N2,  and  produces  the  desired  target  state  I  x0 )  1 10 )  in  0(^lN  )  steps  given  a  randomly  chosen 

telephone  number  to-  We  further  discuss  how  and  why  sequential  (though  not  concurrent  parallel)  searches  can  be 
performed  on  multiple  database  states.  Finally,  we  show  how  this  procedure  can  be  generalized  to  databases  with  more 
than  two  correlated  lists  (e.g.  |x)|t)|s)|r) ...). 


2.  GROVER’S  SEARCH  ALGORITHM 


2.1  The  Grover  iteration  and  the  phase  kickback  or  solution  tagging  operation 

The  essential  operation  in  Grover’s  algorithm  is  the  Grover  iteration’ 

G  =  U  JJ f  (2) 

composed  of  two  functionally  distinct  unitary  components  (i)  U f ,  the  phase  kickback  (PK)  or  “solution  tagging” 
operation,  and  (ii)  U  ± ,  the  inversion  about  the  mean  (IAM).  The  phase  kickback  unitary  operation  works  as  follows 

Uf  | jc) ® | -)  =  Uf  | x) ®  =  | x) ®  =  (_!)/<*)  | ^ 0  I0)"!1)  =  (_!)/<*)  (3) 


where  we  have  used  ( 1 )  with  y)  =  H 1 1)  =  -)  =  (|0)-  l))/>/2  with  the  Hadamard  operator  H  = 


/ V2  which  also 


maps  //|0)  =  |+)  =  (|0)  +  |l))/V2  .  In  the  first  equality  of  (3),  Uf  acts  on  both  terms  of  | yj  simultaneously  (quantum 

parallelism).  Explicitly  working  out  both  casesy(vo)=l  and  f[x^x0)=Q  yields  two  results  that  can  be  encapsulated  into  the 
single  statement  given  by  the  rightmost  expression.  This  is  the  famous  phase  kickback2,3  4  in  which  the  evaluation  of  the 
function/is  stored  in  the  quantum  phases  e,H  (here,  with  9  6  (0, 71}).  Since  the  single  auxiliary  qubit  |  =  |-)  is  returned 

to  its  initial  state  after  the  PK  operation,  one  often  abbreviates  (3)  to 

Uf\x)  =  (-\)Kx)\x)  (4) 
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with  the  understanding  that  the  required  auxiliary  qubit  y)  is  implied.  The  net  result  of  U f  is  that  the  sought  after 
solution  state  |x0)  is  flagged  with  a  -1,  while  all  other  states  |x  ^  xj)  are  unchanged.  An  explicit  matrix  realization4  of  the 
unitary  operator  U ,  is  shown  in  Fig.  1  for  the  case  of  two  qubits  with  solution  state  x0  =  2  in  the  decimal  representation, 

corresponding  to  10  in  the  binary  representation.  (From  now  on  we  will  primarily  use  the  decimal  representation 
x£{0,1,...,N-1}  on  n  qubits,  where  N=2").  Figure  1  illustrates  the  assertion  that  one  needs  to  know  the  solution  x0  in 
order  to  construct  the  oracle  U f  . 


x  =  0  x=l  x=2  x  =  3 
>■  =  0  1  >  =  0  1  >’  =  0  1  >'  =  0  1 


Fig.l  Explicit  construction  of  the  unitary  phase  kickback  operator  UJ  for  the  case  of  two  qubits  labeled  by  the 
decimal  x£{0,l,2,3}  (<->{00,01,10,1 1},  binary)  representation.  In  this  example,  the  solution  state  is  Xq  =  2  (binary, 
10),  X2  denotes  the  2x2  Pauli  crx  bit-flip  matrix  and  I2  denotes  the  2x2  unit  matrix.  The  superscript  xy  on  UJ 

denotes  that  the  PK  operation  acts  upon  the  3-qubit  state  |  x}  ®  |  y)  ,  where  y  is  a  single  qubit  auxiliary  state. 


In  Fig.l,  the  x=  jc0=2  diagonal  block  of  UJ contains  the  2x2  Pauli  bit-flip  matrix  denoted  as  X->  which  flips  the  single 

auxiliary  y  qubit.  All  other  x^x0  diagonal  blocks  UJ  contain  the  2x2  identity  matrix,  denoted  as  I2,  which  leaves  the  y 

qubit  unaltered.  The  superscript  xy  on  UJ  denotes  that  the  PK  operation  acts  upon  the  3 -qubit  states  |  xj  ®  |  yj  .  and  the 

net  effect  is  to  multiply  the  state  |  x0 )  ®  |  -}  by  the  phase  factor  - 1 ,  leaving  all  other  states  unaltered.  Since  there  are 

A=4=22  qubits  in  this  example,  there  are  4  block  diagonals  in  which  to  place  the  bit-flip  operator  X2.  The  choice  of 
which  specific  diagonal  block  is  X2  placed  is  determined  by  the  solution  state  xn.  Formally,  the  PK  operation  has  the  form 

UJ  =  |x0)(jc0 1®  Xy2  ®  II  that  explicitly  illustrates  this  point.  Thus,  the  construction  of  the  PK  operator 

requires  knowledge  of  the  solution  state  x().  This  is  the  primary  reason  why  UJ Ts  given  as  an  “input”  to  the  GSA,  and  is 
considered  as  an  externally  provided  oracle. 


2.2  The  inversion  about  the  mean  operation 

The  second  unitary  in  the  Grover  iterate  (2)  is  the  inversion  about  the  mean  operation,  given  by 

U>f  =H%nU(yH%n  (5) 

where//18”  takes  the  /?-qubit  initial  state|0^to  the  unbiased,  equal  amplitude  product  state  \y/)  =  lj  V/V  ^  'V_n‘  |  x) 
=  //  0).  (the  »-fold  tensor  product  of  H  |0)  of  all  qubits).  For  /7-qubits,  we  will  denote  this  for  simplicity  as 

//|0)  =  |^).  (6) 


The  operator  U  ±  is  defined  as 
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(7) 


UM=\°)> 

Ugl  \x  *  0)  =  -|x  ^  0). 

Note  that  U  does  not  require  knowledge  of  the  solution  state  Xo',  it  simply  flips  the  sign  of  all  states  except  the  standard 
initial  state  x=0.  UQ±  therefore,  has  the  representation  Ugl_  =  |0)(0|-|l}(l|-...-|A-l)(A-l|  =  2|0)(0|-/JV where  we 
have  used  the  completeness  relation  IN  =  ^ 'J_o'|.r)(x| ,  with  IN  being  the  NxN  unit  matrix.  Thus,  using  (6)  we  can 
express  U  _  in  (5)  as 

Uv.  =2\y/)(y/\-IN  .  (8) 

A  straightforward  calculation3,4  reveals  that  U  L  maps  the  amplitudes  cx  of  arbitrary  quantum  state  \(p)  =  ^  cr  |x) 
according  to  cx  I— »  2  c  -  cx  where  c  =2^  cx/N  is  the  average  of  all  the  quantum  amplitudes  cx.  This  is  easily  seen 
since  Avg  =  |^)(^|  is  the  matrix  with  each  entry  taking  the  value  1/A,  that  maps  an  arbitrary  quantum  vector  |$>)to  a 
vector  whose  every  component  is  c  .  Thus,  U ^  performs  an  inversion  of  each  quantum  amplitude  cx  about  its  mean 
value  c. 


2.3  The  Grover  iteration  as  an  effective  rotation  operator 

The  essence  of  the  GSA  is  that  Uf  acts  simply  in  a  2-dimensional  basis  which  we  denote  as  the  2D  column  vector 
(where  T  =  transpose)  \fj)ie,j  >k)..w]  >  and  acts  simply  in  a  different  basis  denoted  as 

|g,.)  =  [| (//),)(//)]  ,  defined  and  motivated  as  follows.  We  can  decompose  the  state  \y/)  (assuming  the  existence  of 

only  a  single  solution  x0  for  simplicity  of  exposition)  as  (using  repeated  index  summation  convention) 

k>  =  VF5|,7+^ll/  =  ^(^5l,4^l'">'cos'k“)+s,n"k'"'>'f''U>;sm','vF'  <9) 

|  y/}  is  defined  as  the  2-dimensional  state  orthogonal  to  |  y/)  defined  as 

|0--sin//|^)  +  cos/?|^). 

The  action  of  Uf  and  U  ±  in  these  bases  are  given  as 


(10) 


Wbad) 

'1  0" 

Wbad) 

Wbad) 

;U  ± 

fk)l 

'1  0 " 

IMl 

r  k)  i 

_K o*)_ 

0  -1_ 

_\v'good}_ 

_~Vgood)_ 

m\ 

0  -1_ 

m\ 

L-k>J 

u . 


and  the  bases  are  related  by  a  simple  rotation  R  in  a  2D  plane 


M 

,k>. 


=  R 


|  V  good  } 

Wm). 


cos  0  sin// 
-sin//  cos// 


tybai) 

I  Wgood 


(11) 


(12) 


From  (11)  we  see  that  in  their  respective  bases,  both  U f  and  U  ±  act  as  the  2x2  Pauli  phase-flip  operator  Z,  =  a. .  The 
application  of  a  single  Grover  iteration  G  can  now  be  carried  out  in  matrix-vector  form  as 

I rm)=<s\¥)=r,ury,\f,)-¥/Jr.  (*")  t/  |g.) 


sV'tGlm\fm) 


(13) 


with 
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(14) 


tC  =W,Gim  =[cos6>,sin6>] 


=  [cos  3(9,  sin  3<9] , 


cos  2  9  sin  2  9 

_-sin2#  cos2<9 

which  rotates  \y/) ,  with  small  initial  angle  0~1  /Vn  from  the  non-solution  \y/bad)  .to  the  new  angle  30  lying 
incrementally  closer  to  the  solution  state  | i//gonil  'j .  After  k  successive  Grover  iterations  |( =  Gk  \  y/) ,  repeated 
applications  of  (14)  yields  y/)k)  =  [cos(2£  +  l)#,  sin(2£  +  l)<9] ,  i.e.  |  ^  =  cos(2^  + 1)<9 1  )  +  sin(2Ar  + 1)6>  |  }  -  The 

GSA  is  terminated  under  the  condition  (for  N»  1)  that  k/2  =  (2k  + 1)9  ~  2k9  ~  2 k/y[N  or  k  ~  I /rVx /4J  (where 


|_xj  denotes  floor(x),  the  nearest  integer  less  than  x).  This  is  the  famous  quadratic  speedup  of  the  GSA  over  the  classical 
brute  force  unstructured  search.  With  these  approximations,  we  have  sin(2k9  +  9)  ~  cos#  =  y/N -l/yfN  and 
cos(2 kd  +  9)  ~  -  sin  9  ~  - 1/ 4n  ,  yielding  |  y/' [k) )  =  -  sin  9  \  y/bad )  +  cos  9 1  y/gooi )  ~  -  l/V/V  |  y/bad )  +  y/N- 1/ y[N  \  y/good ) . 
The  ratio  of  the  probability  to  find  the  good  state  solution  to  the  bad  state  non-solution  upon  a  (von  Neumann) 
measurement  of  the  state  =  G(i)  i//')  is 


P(good)  _  sin2  9 
P(bad)  cos2  9 

In  other  words,  the  initial  amplitude  sinQ  of  finding  the  solution  has  undergone  an  amplitude  amplification  by  a 
factor  of  cot0~V N,  significantly  increasing  the  probably  to  detect  upon  measurement  of  the  final  state. 


yJN-l/yIN 

■[LIn 


■  N . 


(15) 


3.  THE  ADQC  NON-ORACLE  VARIANT  OF  THE  GSA  BY  XU  ET  AL.7 

In  this  section,  we  briefly  review  essential  ingredients  of  the  adiabatic  quantum  computing,  non-oracle  variant  of  the 
GSA  by  Xu  et  al? .  The  motivation  of  these  authors  was  to  circumvent  the  construction  of  Uf  (see  Fig.l)  by  the  agent 
conducting  the  search  (the  searcher)  that  would  require,  as  discussed  in  the  previous  section,  his/her  explicit  knowledge 
of  the  sought  after  solution  state  |x0) .  They  accomplished  this  goal  by  encoding  both  the  names  and  telephone  numbers 

(the  example  Grover  originally  suggests)  into  a  database  state,  which  can  be  constructed  and  stored  at  some  point  in  the 
past,  and  interrogated  later  in  the  future  when  a  random  telephone  number  is  given,  and  the  associated  name  is  sought. 
To  perform  this  database  search,  they  utilized  the  adiabatic  quantum  computing  methodology. 

Briefly,  in  adiabatic  quantum  computing  (AdQC)  the  emphasis  is  on  the  construction  of  a  time  dependent  Hamiltonian 
H(t)  whose  (time-ordered)  exponentiation  yields  a  time  varying  unitary  evolution  U(t).  The  crucial  mathematical  result  is 
that  if  the  evolution  is  slow  enough  (i.e  H(t)/(dH(t)/dt)  “sufficiently”  small  to  avoid  undesirable  level  crossings)  the 
physical  system  remains  in  the  instantaneous  ground  state  i//(l(t)} .  Thus,  one  constructs  a  Hamiltonian  of  the  form  //(/) 

=  1 1  -s(0 1  Hi  +  s(t)  //p.  where  //,  is  the  initial  Hamiltonian  at  time  t= 0  whose  ground  state  is  known  and  easy  to  achieve 
physically,  and  Hv  is  the  problem  Hamiltonian  whose  ground  state  at  time  t=T  is  the  solution  to  the  search  problem.  The 
monotonic  function  s(t)  satisfies  the  simple  requirements  that  .v(O)K)  and  .s(  T)~- 1 ,  as  we  adiabatically  switch  from  //,  to 
Hp  over  the  course  of  the  evolution  t=[0,7], 

Xu  et  al.1  chose  Hl  =  N_*  X,  ,  where  g  is  some  arbitrary  coupling  constant  and  A;  is  the  Pauli  bit- flip  operator  for  the 

;th  qubit.  Since  Xi  |±).  =  ±|±).,  the  ground  state  of  Ht  is  readily  given  by  the  equal  amplitude  product  state 

h(°)>  =  l/ yf~N  ]^[ '  ()'  |  .  The  innovation  of  Xu  et  al. 1  was  to  encode  the  telephone  book  database  into  the  diagonal 

database  operator  D  given  by 

D  =  ZV<I')('I’  (16) 

1=0 
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where  the  name  i  is  encoded  into  the  quantum  state];')  and  the  corresponding  telephone  number  v,  is  encoded  into  an 

(classical)  integer.  The  encoding  of  the  telephone  book  (the  database  consisting  of  a  two  lists,  one  of  names  and  the  other 
of  corresponding  telephone  numbers,  each  of  which  may  be  chosen  as  a  random  permutation  of  the  integers  [0,. .  ,,A-1]) 
takes  place  at  some  time  in  the  past,  and  the  state  D  is  stored  and  saved  until  searching  is  requested.  Upon  the  selection 
of  a  randomly  chosen  telephone  number  v ,  one  constructs  the  problem  Hamiltonian  Hp  as 


Hp  = 


(  jV-1 

5>, 

\i=0 


'X'hVvj  =(^-v,4)2 


(17) 


The  action  of  D  on  an  arbitrary  computational  state  ]/)  is  given  by  £>|z)  =  v,  ]/) .  Consequently,  the  action  of  Hp  on  the 
same  state  /)  is  given  by 


Hp\i}  =  [D-v,IN^  |z)  =  (v,.-v..)  |/)  =  0|/*)=0  iff /  =  /*,  (18) 

indicating  that  the  ground  state  of  Hp  is  given  by  the  state]  /  =  /*) ,  i.e.  the  sought  for  quantum  state  that  was  associated 

with  the  (randomly)  selected  name  v,  in  the  quantum  database  operator  D.  The  rest  of  the  paper  by  Xu  et  al. 7  details  an 

experimental  implementation  of  this  scheme  in  an  NMR  setting.  In  particular,  the  application  of  the  unitary  evolution 
U(t)  is  given  by  usual  Trotter  expansion  of  chopping  up  the  Hamiltonian  II(i)  in  piecewise  constant  terms  Hit )  over 

small  time  intervals  [/,  /+A/]  and  implementing  U (/)  =  U',,(4,A'  which  is  accurate  to  second  order  in  At.  The  key  idea 

of  this  approach  is  to  encode  the  database  into  a  quantum  state,  and  then  subsequently  search  on  one  subspace  of  the 
system  (e.g.  telephone  numbers)  which  is  indexed  to  the  sought  after  item  in  the  other  subspace  of  the  system  (i.e.  the 
corresponding  names). 


4.  GSA  WITH  AN  ENTANGLED  DATABASE  STATE 


4.1  General  amplitude  amplification 

It  is  well  known' 6  that  the  Grover  iteration  (2)  can  be  extended  to  the  more  general  case 

G  =  UvJJf=AU^Uf,  (19) 

where  A  is  the  operator  that  takes  the  standard  initial  state  1 0)  on  n  qubits  to  an  initial  “guess”  state,  often  taken  to  be  the 


equal  amplitude,  unbiased  state  1 1/ /)  =  1  / ^  ^  x) , 


T 

’i/V/v" 

/  \ 

/  \ 

A  0)  =  H 

0 

= 

l/yfN 

=  \v)  =>  A  = 

W) 

0 

_i/V/v 

V  y 

V  / 

(  \ 


\  J 


(20) 


From  left  hand  side  of  (20),  we  see  that  the  initial  state  |0)  picks  out  the  first  column  of  A,  which  is  |  y/)  .  The  rest  of  the 
columns  of  A  can  be  chosen  arbitrarily,  only  under  the  restriction  that  A  is  unitary  A4I=/N,  requiring  that  all  columns 
(and  all  rows)  are  mutually  orthonormal.  Note  that  in  standard  Grover  iteration  (2)  A  =  H  on .  Further, 
A 1 0)  =  1 1//)  ensures  that  U  ±  =  AU^A1  =  IN  is  again  the  inversion  about  the  mean  operator  (8). 


4.2  The  quantum  database  state  and  the  subspace  phase  kickback  operation 

We  now  develop  a  quantum  database  state  \y/jb)  ,  this  time  in  the  quantum  circuit  model  approach  in  which  we  explicitly 

state  the  unitary  evolution  operators  (vs  the  AdQC  approach,  in  which  the  focus  is  on  the  constructed  Hamiltonians). 
Again  we  utilize  the  example  of  a  telephone  directory  database.  We  will  encode  both  the  names  and  the  telephone 
numbers  into  quantum  states,  and  illustrate  our  implementation  explicitly  with  the  example  utilizing  n= 2  qubits  for  the 
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names  and  n= 2  qubits  for  the  telephone  numbers,  while  concurrently  developing  formulas  for  an  arbitrary  number  n  of 
qubits.  We  consider  the  case  of  N=  2"  (name,  telephone  number)  pairs  .  The  quantum  database  state 

k,»)is  given  by 


(21) 


which  is,  in  general,  an  entangled  state  between  the  name  and  telephone  component  states. 

Note  that  \y/dh)  is  an  Wdimensional  vector  in  a  N2  dimensional  Hilbert  space,  where  the  most  general  state  is  given  by 


1  W-1JV-1  1  jV2-l 

'*')  -  -7=fZI«iy  ®  I A  s-7=t2>*  \k)x,  e  H*  2 


yfN*f 


■0  j= 0 


(22) 


whereHt  and/-/(  are  the  /V-dimcnsional  Hilbert  spaces  of  the  names  and  telephone  numbers,  respectively.  In  (21)  and 


(22)  we  use  a  subscript  notation  to  denote  which  Hilbert  space  the  ket  belongs  |i)  eH(,  |  j}  e  H  t  and  k)a  e  H  X®H  , . 
Let  us  consider  the  specific  example  of  n= 2,  N=  2"  =4,  and  utilize  the  decimal  representation  of  the  states  (i.e. 
{X',ti}ie[ o  i2  3])-  Consider  a  telephone  directory  and  corresponding  database  state  given  by 


names  telephone  #s 
0  2 

{Xi’/'},s[  0,1, 2,3]  =  *  ^ 

2  0 

3  1 


l^>=^(l°)J4+l1)J3)(+l2)J0)(+l3)J1),) 


(23) 


=^(I2L+I7L+I8L+I13U» where  \xi)M),^\Nxi+ti)*r 


Note  that  while  we  have  ordered  the  names  in  (23)  sequentially,  the  two  lists  of  names  and  telephone  numbers  can  in 
general  be  chosen  as  random  permutations  of  the  integers  [0,1,..., Wl].  Our  rational  for  constructing  the  database  state 
\ysdh)  is  simple.  Given  the  telephone  directory  (the  database)  we,  as  the  eventual  searcher  (database  interrogator),  can 

encode  this  classical  information  into  the  quantum  state  |  y/dh )  and  store  it  for  subsequent  interrogation.  Suppose  at  a  later 
time,  we  select  (or  are  provided  with)  a  random  telephone  number  t,  and  desire  to  find  the  associated  corresponding 
name  x*.  We  can  then  construct  the  phase  kickback,  telephone  number  tagging  operator  U\y  utilizing  the  known 


information  of  the  selected  telephone  number  t.  For  example,  if  t=  2,  the  operator  U\y  would  have  the  form  given  in 
Fig.l  (with  x  now  replaced  by  t ).  Note  that  U‘f  acts  on  the  ty-subspace  (indicated  by  the  superscript)  of  telephone 


numbers  t  and  the  auxiliary  qubit  y,  and  not  on  the  x-subspacc  of  names,  on  which  we  are  seeking  the  associated  name 
x*.  Thus,  in  the  full  2A^-dimensional  Hilbert  space  of  H  r(W)  <8>  H  t(N)  ®H  v(1)  (where  the  superscript  denotes  the  dimension 


of  the  Hilbert  space)  the  telephone  number  tagging  operator  has  the  following  form,  and  operational  PK  effect 

\l  t  =  t 
I  0  t  &  t 


Ux»=rN®U‘fl  where  ft.(t)  =  8.= 


rN®u'i\x)®\t) 


|0),-|1) 

V2 


I0).-!1). 

^  ' 


V2 


n  (24) 
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Note  that  U'/t  performs  an  effective  sign  flip  on  states  |x)v  ®  for  all  values  of  x.  Due  to  the  tensor  product  nature  of 
the  component  states,  a  PK  sign  flip  on  |  f'j  produces  an  effective  PK  sign  flip  on  |x^  ®  |t*)  ,  which  includes  the  sought 
after  state  |  x  ^  ®  1 f  ^  .We  now  proceed  to  construct  the  Grover  iteration  in  this  approach. 


4.3  Encoding  the  database  into  the  quantum  database  state 


From  (20)  we  need  to  construct  a  unitary  operator  A  =  Axl  such  that  ^4 1 0^  ®  |0)  =  4|0)  t  =  | y/db . 

4°>,®l°)i  (25) 

This  is  most  easily  accomplished  if  we  perform  a  relabeling  of  the  indices  of  the  xt  component  states  in  the  lower, 
rightmost  line  of  (23)  so  that  we  bring  them  to  the  first  N  entries  of  the  N2  database  vector,  i.e. 

k,ft)-l/V4(|2)„+|7)H+|8)H+|l3)ti)^  1/V4(|°'L+l1'L+l2'L+|3'U^k;6),  which  we  will  call  the  prime 

frame,  which  we  denote  by  primes  on  the  component  values.  In  the  prime  frame,  the  database  state  has  the  form 
\y/'db)  =  l/  s/4  [1,1,1, 1,0,...,  0]r  and  we  seek  a  unitary  operator  A'  with  the  property  that  4'|0')  (  =  |  y/'dh  .  Though  there  is 
much  freedom  in  choosing  such  an  A! ,  the  simplest,  most  direct  (though  non-unique)  choice  that  we  adopt  here,  is  to 
choose  A'  ~  Hn  =  II" ,  the  «-fold  tensor  product  of  2x2  single  quibt  Hadamard  unitaries.  In  the  prime  frame  A'  takes 
the  block  diagonal  direct  sum  form, 

A'  =  Hn®In1_n,  (26) 

(see  Fig.  2),  in  which  7^  is  the  ( N2 -N)  x  (N2-N)  identity  matrix  acting  on  those  states  |i)  ®  |i),  not  in  \y/'db)xl  ■ 


N 


N 


N*-N 


Fig.2  Form  of  the  unitary  A'  operator,  effecting  the  operation  A'  1 0f}r(  =  |  y/'db ) ^  ,  in  the  prime  index  ordering  for 

which  the  N2 -dimensional  primed  database  vector  |  Wdb)xt  has  equally  weighted  values  in  the  first  N components. 

Hn  is  the  n-fold  tensor  product  single  quibt  2x2  Fladamard  unitaries.  Ia^-jv  is  the  (N2-N)  x  (N2 -N)  identity  matrix 
operating  on  those  states  1 1)  <S> | _/)  not  in  \y/'db)  t  ■  A1  has,  the  block-diagonal  direct  sum  form 

a'  =  hn®in2_n. 

We  can  transform  A'  back  to  the  original  unprimed  frame  by  a  series  of  N2  x  N2  unitary  operations  St  .  =  .S',1'  that  swaps 
rows  i  and  j  of  any  matrix.  In  our  particular  example  (23)  where  |  y/db)  =  l/x/4  (|2)n  +  \l)xt  +|8)w  +|l3)vJ  we  have 

k*)  =  ^(|2)»  +  |7)a  +|8)I(  +|l3)i()  =^=Y^k^xi\  {k^k^k-,,^)  =  {2,7,8,13} 

k*)=^(|°'L +I1'),* +I2'L +I3'L)-^EI''L’ 

-Tls  A'-S  S  S  S  A'  (S  )  -I1  0’0e{0'0),0,0} 

U Sf*  -  °'-2  r'7  2'-8  3'13  ’  W ~  0  otherwise 


(27) 


Proc.  of  SPIE  Vol.  8057  80570R-8 


Downloaded  from  SPIE  Digital  Library  on  21  Feb  2012  to  132.3.29.68.  Terms  of  Use:  http://spiedl.org/terms 


which  is  illustrated  in  Fig.  3.  The  swap  operators  Sy  acting  on  a  quantum  state  vector  effectively  performs  a  Pauli  bit- flip 
operation  X2  between  the  ;th  and  /ill  components,  and  therefore,  acting  on  a  matrix,  Sy  swaps  the  /ill  and  /th  rows.  Note 
that  we  perform  the  row  swaps  from  the  prime  to  the  unprimed  frame  beginning  with  the  largest  value  of/'  =  7V, 
backwards  to  smallest  value  i'  =  0 . 


4.4  Construction  of  the  Grover  iteration 

The  construction  of  the  Grover  iteration  is  most  clearly  described  in  the  prime  frame.  However,  in  the  numerical 
examples  discussed  below,  the  simulations  are  carried  out  in  the  unprimed  frame.  As  a  (common)  slight  abuse  of 
notation,  we  will  drop  the  explicit  labeling  of  the  auxiliary  y  qubit  (leaving  it  implied)  that  is  necessary  to  perform  the 
PK  operation  in  the  telephone  number  subspace.  We  will  denote  this  unitary  operator  simply  as  U\  enacting  the 

effective  conditional  sign  flip  TJ‘f%  |t}(  =  (— analogous  to  (4).  We  shall  further  write  the  full 


3'  4)  «?  !9#  11C  D  B  15 


VA  U  l/A  1 /A 
i/A  -i/4  l/A  A4 
i/A  u  -VA  -¥A 
i/A  -I/A  -I/A  m 


0 


A' 


0 


=> 


0 
1 
2 

3 

4 

5 

6 

7 

8 

9 

10 
11 
12 

13 1/4  -i/4  -i/4  i/4 

14 

.V-l 


i/4  i/4  i/4  i/4 


i/4  -i/4  i/4  A4 
i/4  V4  -i/4  -V4 


456  789  10  1112131415 
1  ] 
% 


A:] 


^  -  ]”[  —  S0,2SY1S2,8S3l3A 


i= 0 


Fig.3:  Successive  row  swapping  operations  to  transform  A'  in  the  prime  frame  to  A  in  the  unprimed  frame  for  the 
specific  telephone  database  example  in  (27). 


N2  x  N2  operator  in  (24)  acting  on  xt  vectors  as  Uxl  =  IXN  ®U'r  .  This  matrix  is  block-diagonal  in  form,  which  each  block 
consisting  of  the  NxN  sub-matrix  telephone  number  tagging  operator  U\  .  Formally  Uf  takes  the  form 


t/"  =  rN  ®  u‘fm  =  ®  u‘f'  =  u‘u 


/,*  ’ 


(28) 


which  is  the  tensor  product  of  the  NxN  matrix  U'f  acting  in  the  TV-dimensional  /-subspace,  for  a  given  value  of  x  (for 


which  there  N  possibilities),  with  the  identity  matrix  /*  in  the  x-subspace.  The  direct  product  form  of  (28),  and  presence 
of/^j  in  the  x-subspace  indicates  that U‘f_  transforms  |x^  ®|/*^  >-|x^  ®|/*^  (leaving  all  other  states  unchanged), 

unbiasedly  for  all  value  ofx,  not  just  for  the  sought  for  value  x=x  . 


Turning  to  the  N2  x  N2  inversion  about  the  mean  (LAM)  operator  U"  we  note  that  both  the  states  |  O') u  and|0)n  in  the 
prime  and  unprimed  frame,  respectively,  take  the  form  of  the  vector  [l,0,...,0]r .  Therefore,  the  operator  U'*1  has  the 
analogous  operational  effect  as  in  (7) 
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where  use  has  been  made  of  the  completeness  relation  /  2  =  .  (k  | ,  with  I,,2  the  N2xN2  unit  matrix.  Further, 

in  the  unprimed  frame  U  "  has  the  exact  same  form  as  U"‘  without  the  primes.  In  both  the  primed  and  unprimed  frame 

0X  a1 

we  can  construct  U'xt  or  U*  independent  of  knowledge  of  t  and  x  . 

o-1-  ox 

Using  the  fact  that  A'  |0')tf  =  |  ¥'db)xt  yields 

=2k*L  xl(V/db\~^N2, 

s(2KL>^-4)©(-^J,-  (30) 

=  u^©((-4)©...©(-4)). 

V  v - V - ' 

N-l 

It  is  important  to  note  that  |  y/'db)xl  (y'db  |  is  the  averaging  operator  (see  the  discussion  following  (8))  only  on  the  first  N 
components  of  a  primed  quantum  state  vector,  with  the  effect  of  averaging  only  the  components  of  the  database  state. 
Even  in  the  unprimed  frame  where  the  N  components  of  \ysdb)xl  arc  spread  out  over  the  N2  possible  positions, 

|  \(/dh  )w  i7  (y/dh  |  sti  1 1  performs  an  average  only  on  the  N  components  of  \y/db)xt-  Thus  the  N  x  N  matrix 
UT  s  2  K )„  x,  (v*  \-  /v  performs  an  inversion  about  the  mean  operation  on  the  first  N  components  of  an  N2 
component  primed  xt  quantum  state  vector.  (The  superscript  on  U'‘ ' 1  denotes  the  dimension  of  the  subspace  on  which  it 

acts).  The  full  Grover  iteration  (19)  is  then  given  by 

Gxt  =  Ux\  Uxt. , 

v  f 

=  Ax‘  U“  Axr'Vxl. ,  (31) 

=  (2k<ft}_,K|-40(/«®c/>)’ 

which  we  illustrate  in  Fig.  4  in  the  prime  frame  (for  clarity)  acting  on  |(/4>)„  • 


.v 


N 


N 
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Fig.4:  Illustration  of  the  action  G'xt  =  U’x[U'xt.  =  (^-\y/db)xt  „  {vdb  \  _42  )|^iv  ®  T4  )  >  the  Grover  iteration 
(31)  in  the  primed  frame,  acting  on  the  |  \j/'dh )  f ,  the  primed  database  state.  Note  that  in  each  NxN 
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block  t/L  performs  a  sign  flip  U' .  \t)t  =  (—  l)^  conditioned  on  the  selected  telephone  number  t*, 

independent  of  the  value  of  x.  The  operator  U'i±)  =  2|  w'db)x,  „  (Vm  \  ~ Tv  (where  we  have  suppressed  the 
subscripts  in  the  figure)  performs  an  inversion  about  the  mean  on  the  N  components  of  I  y/[1b )  . 


The  net  effect  of  (31)  is  that  on  the  first  N  components  of  the  primed  database  state  |  y/’db )  _ ,  we  affect  a  Grover  iteration 
of  the  original  form  in  (2).  On  the  later  N2-N  components  of  |  ¥db ) , we  perform  the 
operation  |  x  ^  x  'j  ®  1 1 *  ^  - 1  x  x  'j  ®  | ?* )  in  each  NxN  block,  followed  by  the  multiplication  by  the  NxN  identity 

matrix  Iw.  However,  since  the  latter  N2 -N  components  of  |^6)  are  initially  zero  (which  we  will  generalize  below),  they 
remain  zero  after  the  Grover  iteration.  Thus,  after  k  ~\tT'Jn  /■  4  J  Grover  iterations,  the  amplitude  of  the  state 
|x*^  ®|^*}  tying  somewhere  in  the  first  N  (of  the  N2)  components  of  y/’dh the  will  be  driven  to  a  magnitude 
-s/iV  —  l/ ~Jn  ~  0(1)  (see  the  discussion  following  (14))  with  probability  0(1  -  l/N)  for  detection  upon  measurement. 


4.5  Numerical  simulations  of  algorithm 

In  Fig.  5  we  show  a  simulation  for  the  case  of  n= 3  qubits  where  a  pair  of  arrays  of  names  and  telephones  of  size  N=2"=8 
are  chosen  as  random  permutations  of  [0,1,. . .,7V-1].  In  the  code,  the  database  state  | y/db)r,  n'|^)w  is 

constructed  in  the  unprimed  frame,  and  from  the  specific  collection  of  jV  indices  { k } ,ih  in  the  database  state,  we  construct 
A  in  (27)  from  A' ,  as  discussed  in  (26)  and  illustrated  in  Fig.  3.  This  allows  us  to  construct  U‘\  .  We  next  generate  a 

random  telephone  number  t  and  use  it  to  construct  the  specific  PK  telephone  tagging  operator  U'[.  .  Subsequently,  we 

assemble  the  Grover  iteration  Gxt  =  t/”  Uxt.  and  apply  it  for  \jt4n /4J  =  2  .  In  Fig.  5  note  that  only  N=  8  of  the  total 

;V2-64  probabilities  are  non-zero  throughout  the  whole  evolution,  corresponding  to  the  N=  8  non-zero  amplitudes  of  the 
database  state  | y/db)xt  ■ 
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iteration  =  1  of  2,  prob(max)  =  0.7812 


iteration  =  2  of  2.  prob(max)  =  0.9453 


Fig. 5:  Grover  iterations  (31)  for  n= 3  qubits  (iV=23=8)  with  a  randomly  selected  telephone  number  t  =  2.  The  plots 
(left  to  right)  show  the  probabilities  (amplitude  squared)  for  the  {0,1,  |^;rV/V  /4 J  =  2}  iterations.  The  abscissa  is 

the  combined  xt  index  k=Nx+t  ranging  from  0  to  A2-l=63.  The  Grover  iterations  acts  on  the  | )  and  drives  it 

towards  the  state |/c  =  10^  =  |x  =  1^  ®|t  =2^  with  near  unit  probability.  Note  that  only  N=  8  of  the  total 

A2=64  probabilities  are  non-zero  throughout  the  whole  evolution,  corresponding  to  the  N=  8  non-zero  amplitudes 
of  the  database  state  |  y/db  )  . 


5.  CONSIDERATIONS  AND  EXTENSIONS 

5.1  Initial  state 

It  is  illuminating  to  consider  the  action  of  our  Grover  iteration  Gxt  on  initial  states  |  y/lml )  other  than  the  constructed 
database  state  |  Ydb)xl  •  In  general,  the  normalized  initial  state  could  be  written  in  the  form 

I ¥inU )xt=\[p  I Vdb )xt  +  yfr-P  I ¥ndb )xt  >  (32) 

where  |  y/n4b )  _  denotes  a  normalized  non-database  state,  i.e.  the  state  formed  by  all  components  not  in  the  database  state 

|  ¥db)xt  ■  3*1  (32),  p  =  |  (ydb  |  ¥  mu  )P  *s  the  probability  to  find  the  initial  state  in  the  database  state.  After  the  |^;rV/V  /  4  J  of 
Gxt  the  final  state  is  approximately 

1^)*  ~  yfp\x}y)l  +  'fi-p\¥**)*’  (33) 

which  implies  there  is  only  a  probability  p  to  detect  the  sought  after  solution  state  |  x*  'j  1 1“  'j  .  Thus,  as  long  as  p>l/2,  the 
form  the  GSA  presented  here  does  better  than  its  classical  0(N/2)  exhaustive  search. 

The  initial  state  (32)  might  occur  as  an  imperfect  attempt  to  construct  the  desired  database  state  | ¥db)x,  ■  A  simpler  state 
to  form,  is  the  N2  equal  amplitude  state  \¥ni)  =#f”  ® |0}r  ® |0),  since  the  last  equality 

shows  that  we  can  form  this  state  directly  by  the  application  of  2«-fold  tensor  product  of  Hadamards  H ®”  ®  Hfn  acting 
on  the  tensor  product  of  the  n-qubit  standard  name-state  1 0)  and  the  n-qubit  standard  telephone-state  |0)  .  However, 
from  (32)  and  (33)  such  an  initial  state  renders  the  GSA  worse  than  classical  exhaustive  search  since  there  are  N2 -N non¬ 
database  states  each  with  probability  1/ N2  that  are  unchanged  by  the  Grover  iteration  Gxt  for  a  total  probability  for  the 
detection  of  \ip„db)xi  upon  measurement  given  by  1-1/A.  The  initial  probability  of  p=\!N  to  find  \¥N2/  i11  the 
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state  | Wib)xt  remains  the  final  probability  to  find  y/flnal 'j  in  the  solution  state  |.r*^  .  Figure  6  illustrates  though  that 

G'1  the  does  act  only  on  the  )V-component  state  y/dh )i(  buried  within  N 2  sized  initial  state  \yN2  )  • 
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Fig. 6:  Grover  iteration  (31)  for  n= 3  qubits  (7V=23=8)  with  a  randomly  selected  telephone  number  t  =• 6  with  initial 
state  | =  \  ys Ni  'j  =  l/ ' \lN2  y„.  \k)xt  ■  The  plots  (left  to  right)  show  the  probabilities  (amplitude  squared) 

for  the  {0, 1,  |^;rV/V / 4 J  =  2}  iterations.  The  abscissa  is  the  combined  xt  index  k=Nx+t  ranging  from  0  to  N2- 1=63. 
The  Grover  iterations  acts  upon  the  \vdb)xt  portion  of  |  Vinit  )xl=y[p  \  Vib  )x,+'J1~p\ Vmb ),, where 

p  =  l/yjN2  =1 /N  ,  and  drives  it  towards  the  state |&*  =  30^  =  |x*  =3^  ®|t*  =  6j  with  probability  p=l/8. 
Note  that  all  7V2=64  probabilities  are  non-zero  throughout  the  whole  evolution,  but  only  the  N=  8  amplitudes  of  the 
database  state  Wdb),,  are  acted  upon  by  Gxl  (compare  with  Fig.  5).  Because  the  amplitudes  in | )  are 

unchanged  by  Gxt  using  the  initial  state,  |  y/inH  )  =  |  y/Ni  )  yields  inferior  performance  when  compared  to  classical 
exhaustive  search. 


5.2  Sequential  vs  concurrent  searches  with  a  generalized  database  state 

Since  the  essence  of  the  above  variation  of  the  GSA  is  to  implement  the  Grover  iterations  in  an  (V-dimensional  subspace 
of  a  larger  A2  dimensional  Hilbert  space  (of  xf-states)  it  is  natural  to  inquire  if  multiple  independent  Grover  search  could 
be  executed  in  M<N  independent  subspaces  concurrently.  The  answer  is  unfortunately  no,  for  a  reason  illuminated  by 

the  form  the  our  current  algorithm.  Consider  |  'V  DB )  =  i/a/tT  ^  ^  1 'j  =  l/^M  N  ^ Q'  ^  ^  |  x\m)  'j  ®  |  t\m)  ^  an  M 
database  state  encoded  containing  MN  <  N2  elements.  We  assume  that  a  given  pair(xl('"),tl(m)) ,  the  state  ®  |  t\m) 'j 

appears  in  one  and  only  one  database  state  |  t//d" 1  ^  .  The  crucial  operation  is  the  t-subspace  PK  operation  Uxf[  (28), 
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which  is  applied  independent  of  the  x-subspace  Ux!  =  IXN®U\  =  ®  U',  =  U‘,  .  Suppose  M=  2,  and  we 


attempted  to  apply  two  simultaneous  ,Vx,V  PK  operations  U'f,  and  U' .  given  two  separate  telephone  numbers  f°r 

Jtx  Jt2 

two  different  database  states  Wdb)  and  ■  This  would  imply  that  we  would  have  to  implement  an  operator  of  the 

form  Uxt. r,  =  U[.  ©  U[,  ©C/ll2  „„  (where  in  the  prime  frame,  f/'  ,  acts  on  the  first  N  components  of  consisting  of 


kk 


k 

N 


k 

N 


N'-IN 

N2-2N 


Wdb)  >  U‘.  acts  on  the  second  N  components  of  I  xf/DS)  consisting  of  Wdb)  and  4  ,  is  some  arbitrary  unitary 
operator  acting  on  the  remaining  zero  valued  N-N  components  of  |'f/Dg))-  However,  the  construction  of  such  a  PK 
operation  requires  some  specific  knowledge  of  the  x-subspace  (i.e.  about  location  of  sought  for  solutions  {x*,x*}  ),  and 
therefore  would  not  of  the  x-egalitarian  tensor  product  form  IXN  ®  U‘  .  ,  ,  where  now  U' ,  .  acts  solely  in  the  ?-subspace, 
independent  of  the  values  of  x. 

The  best  approach  is  to  apply  IXN®U‘.  and  its  appropriate  averaging  operator  U1'  =  \2\i//(l!)  (1//$  \-I  2  (until  we 

Jtx  w.  \  I  /  Xt  Xt  \  I  N  ) 


converge  Wdb 


to  the  solution 


t*)  ,  and  then  repeat  the  process  applying  /*  ®  U‘ .  and  its  appropriate 


.  The  case  of 


averaging  operator  Ux\  =(  2 \y/™ )  ly/jf  -/  ,  )  until  we  converge  Wdb)  to  the  solution  xj)  ®  M 

y/ j  \  1  I  Xt  Xt  \  I  ty  ]  I  I  Xt  I  /  X  I  “  /  t 

M=  2  is  illustrated  in  Fig.  7,  in  the  prime  frame,  where  y/f )  occupies 
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Fig. 7:  Illustration  of  the  sequential  GSA  in  M=  2  separate  N-dimensional  subspaces,  illustrated  in  the  prime  frame. 
The  two  N-dimensional  subspaces  Wdb)  and  Wdb)  occupy  the  first  and  second  set  of  N  components  of  the 
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joint  database  state  I'F^)  —  i^V^db  )  t  +  |^/'(*))  / ) /-\/2  .  (a)  Given  two  telephone  number  (tx  ,t2 )  for  the  two 
databases  Wdb)  and  ,  the  operations  Ux{  =rN®U'f  and  Ux'±  =(2\l//^\  {¥db\~ I (the 

same  in  Fig.  4)  drive  |  A1  DB )  — >  (\xl  \  0 1  t,  \  + 1  \  )  42  .  (b)  Subsequently,  constructing  C/j/  =  I*N  0  U'ft 

\>  '  X  \  '  t  \  I  Xt  /  /  12  12 

and  using  Ux\  =  (ipff  |  -/^  j  drives  the  solution  to  — »  (j^*)  .®|^i*)  +  |*2)  ®\4)^ j 42  . 


components  0  to  At- 1  of  |xl/DS},and  \¥df)  occupies  components  At  through  2N-1.  Since  both  \y/jb)  and  |  Vm)  have 
been  encoded  sometime  in  the  past  into  separate  N-dimensional  subspaces,  we  know  operators  A[  and  A[  that  encode 
| ¥db)  hd0  the  first  N  components  of  in  the  prime  frame,  and  encode  |(^',*))  into  the  second  N  components. 

Hence  we  know  how  to  construct  the  averaging  operators  Ux\  =  A[U\AX'  and  Ux\  =  A^U'.A £  that  act  on  the  first  and 
second  set  of  N  components  respectively  in  I'F^)  •  The  net  effect  of  the  sequential  GSA  as  illustrated  in  Fig.  7  is  to 
drive  to  the  joint  solution  =  (l^*)  0  |0jyV2  ,  with  probability  1/M  to  measure  a  specific 

solution  xm  ®|/*J  ,  m£  { 1,2}.  Statistically,  this  is  equivalent  to  performing  two  separate  GSA  on  two  separate 
database  states  | ¥db)  and  |^>)  as  originally  discussed  in  Section  4.4  above  and  illustrated  in  Fig.  4.  The  only 

advantage  of  performing  sequential  GSA  as  illustrated  in  Fig.  7  is  to  continue  using  the  original  database  state  |VFDS^, 
without  re-creating  it  for  subsequent  subspace  searches. 


5.3  Formation  of  the  initial  state 

Once  the  telephone  book  database  is  given,  the  quantum  database  state  is  constructed  via  (25)  from  the  unitary  operator 
A  and  the  standard  initial  state  |0)  ®|0)  =|0)  byH|0)o  =  | y/  db )vj .  In  this  section  we  briefly  show  that  the  initial 


N- 1 

database  state  =  ®|tf)  given  in  (21)  can  also  be  constructed  in  0(log  N)  =  O(n)  cascaded  stages 

i=0 

(using  N=2"  for  simplicity),  where  parallel  operations  can  be  utilized  at  each  stage  This  assertion  rests  on  the  following 
observation.  Consider  two  arbitrary  states  |  </>)  and  |  (p ) .  Define  the  states  I ±)  =  (U)  ±  I  <pj)  j 42  and  the  Hadamard  operator 
defined  on  these  two  states  as  H.  =  |  +)  (tfi  |  + 1  -)  (ip  |  with  the  effect  |  =  |  +)  and  H ^  |  <p}  =  |  -) .  A  simple  calculation 

shows  that =  [\<!>) (<!>\  +  \(p) (cp\)  +  {<j)\(p)\+) (-\  +  (^)|-}(+|,  which  is  the  identity  operator 

If  =  +  |^)(^|  °n  the  subspace  of  these  two  states  if  and  only  if  |^)and  |#>)are  orthogonal,  i.e.  =  In 


other  words,  II 0  takes  the  form  of  the  usual  Hadamard  operator  —= 

4  2 


1  1 
1  -1 


if  the  rows  and  columns  are  labeled  by  the 
®|t,.)(  in  |^rfA)are  orthogonal  (and  known 


orthogonal  states  |^)and  |  (p)  .  Since  each  of  the  component  x-t  states  xt  ®  |  /■ ) 
since  we  are  initially  encoding  the  telephone  directory  database),  we  can  construct  the  necessary  Hp  to  combine  them 
pair  wise  in  a  logarithmic  number  of  stages  to  successively  construct  |  y/db )  .  This  is  illustrated  in  Fig.  8  for  the  database 
state  K)  =  l/V4(| 0)x  1 2\  +  \ l)x  | 3),  + 1 2\  | 0),  + 1 3)_r  1 l), )  given  in  (23). 
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|0)J2>, 

I1),!3), 

l2>J°>, 


Fig. 8:  Construction  of  database  state  \Wdb)  =  l/>/4  (|0^  |2^(  +  |l^  |3^(  + 1 2^ |0^  + 1 3^  |l^j  by  successive 

Hadamard-type  operations  on  successive  pairs  of  orthogonal  states  (see  text)  for  the  case  of  n= 2  qubits,  i.e.  two 
N=22=4  dimensional  subspaces  to  represent  the  name  (x-states)  and  telephone  number  (t-states)  subspaces 
respectively.  For  general,  N=  2"  dimensional  x-  and  t-subspaces,  the  first  (leftmost)  stage  has  N/2  gates  Hxt  acting 

on  known  (since  one  is  encoding  the  telephone  directory  database)  pairs  of  orthogonal  input  states  ®  j}  and 

|  i'}  ®  |  j'}  acting  in  parallel.  At  each  subsequent  stage  (moving  rightwards)  half  the  number  of  Hadamard  gates 
are  constructed  and  used  to  form  further  superpositions  of  orthogonal  input  states,  which  successively 
constructs  |  y/db  )  after  0(log  iV)=0(n)  stages. 

6.  SUMMARY  AND  CONCLUSIONS 


The  focus  of  this  paper  was  to  illustrate  a  variant  of  Grover’s  algorithm  for  the  case  of  a  k=2  indexed  database  state  of 

N- 1 

the  form  | y/db)  ,  =  1  /\N  ’yjjt,-)  ®| t,)d .  The  rationale  our  approach  was  to  avoid  the  requirement  that  the  oracle, 

(=0 

implementing  the  phase  kickback  operation,  has  to  be  supplied  to  the  searcher  of  the  database  by  an  external  agent  (as  in 
the  original  formulation  of  the  GSA).  Instead,  our  variant  of  the  GSA  was  designed  so  that  the  searcher  could  initially 
encode  the  database  state  and  subsequently  search  it  at  a  later  time,  without  having  to  know  the  sought-for  result  in  order 
to  construct  the  phase  kickback  operator. 


The  variant  of  the  GSA  discussed  in  this  work  can  easily  be  extended  to  multi-indexed  databases  states  of  the  form 

N- 1 

Wdb)mr  =  l/V/V ^|x,.)  ®|t(.)  ®|s,)t  If  this  generalized  database  state  is  encoded  in  the  past,  then  a  later 

i= 0 

time  a  chosen  subspace  component  (e.g.  the  /-subspacc  telephone  number  as  illustrated  in  this  paper)  can  be  searched  on 
through  the  construction  of  a  phase  kickback  operation  U'f  =  IXN  ®  U'f  ®  7^  ®  I'N  on  that  subspace,  implementing 

f[t*)=  1  andy(?^t*)=0  for  a  given  t,  producing  the  result  U'f  |x,.)v  ®  |t,.)(  ®  |j()  ®  |^.)(  =  | jcf ®  (_K)<)  ® \s‘)s  ®  | ri)r- 
=  -|xf)r  ®  |t,-)(  ®  |s,.)(  ®|^)(.  For  a  ^-component  database  state  (i.e.  k  different  index  states  |jq)  ,|/,)  ,|sf)  , | /j )  ,...  ), 
general  amplitude  amplification  as  given  in  (19)  G  =  U ^  U f  =  AUl}l  A  'Uf  can  be  used  to  perform  an  0(\[N  )  Grover 

search  algorithm  in  the  N  dimensional  subspace  of  a  general  (V4  dimensional  Hilbert  space.  Again,  A  is  the  unitary 
operator  that  takes  the  standard  state  |0)  to  the  database  state  |^r)  .  The  utility  of  this  approach  depends  upon  the 

ease  and  efficiency  of  constructing  the  operator  A.  and  hence  the  quantum  database  state  (//)  ;  .  Similar  work  along 

the  lines  of  a  quantum  Grover  search  upon  multi-index  states  has  been  considered  by  Pang  et  al.s. 


DISCLAIMER:  Any  opinions,  findings  and  conclusions  or  recommendations  expressed  in  this  material  are  those  of  the 
author) s)  and  do  not  necessarily  reflect  the  views  of  AFRL. 
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